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Abstract 

Electron spin resonance studies of the Ni-doped spin-Peierls compound CuGe03 has been per- 
formed in the frequency range 9-75 GHz and temperature interval 1.3-20 K. An anomalous tem- 
perature dependence of the ^-factor below the spin-Peierls temperature was observed for doped 
samples. At low temperatures the g-factor is much smaller than the value expected for Cu 2+ and 
Ni 2+ ions and is much more anisotropic than for an undoped crystal. This anomaly is explained 
by the formation of magnetic clusters around the Ni 2+ ions within a nonmagnetic spin-Peierls 
matrix. The formation of magnetic clusters is confirmed by the observation of a nonlinear static 
magnetic susceptibility at low temperatures. 

The reduction of the spin-Peierls transition temperature was found to be linear in the dopant 
concentration x in the range < x < 3.2%. The transition into the antiferromagnetically ordered 
state, detected earlier by neutron scattering for x > 1.7%, was studied by means of ESR. For 
x=3.2% a gap in the magnetic resonance spectrum is found below the Neel temperature and the 
spectrum is well described by the theory of antiferromagnetic resonance based on the molecular 
field approximation. For x=1.7% the spectrum below the Neel point remained gapless. The 
gapless spectrum of the antiferromagnetic state in weakly doped samples is attributed to the small 
value of the Neel order parameter and to the magnetically disordered spin-Peierls background. 
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1 INTRODUCTION 

The magnetic properties of crystals of the quasi-onedimensional magnet CuGeOs have 
been extensively studied since Hase et al JTJ reported that this compound is the first 
inorganic spin-Peierls material. The spin-Peierls transition occurs because the S = 1/2 
Heisenberg antiferromagnetic chains are unstable when coupled to a three-dimensional 
phonon field [0. Below the transition temperature the magnetic chains are dimerised 
and the distance between the neighbouring magnetic ions as well as the exchange integral 
alternate. The positions of pairs of magnetic ions after the dimerisation are correlated 
between neighbouring chains. Thus the dimers construct an ordered sublattice. 

The dimerization of the magnetic ions in chains arranged along the c-direction of the 
orthorhombic crystal leads to the formation of the non-magnetic ground state separated 
from the excited triplet states by 2 meV ~ 23 K [[| . The unit cell in the 

dimerised state is doubled along a and c directions and the intrachain exchange integral 
takes the alternating values J\p = (1 ±5)J C . Here J c =10.2 meV is the intrachain exchange 
integral and 5 ~ 0.04 is a distortion parameter ||]. The existence of such dimerization was 
confirmed by the observation of additional reflections by X-ray || and neutron Q diffrac- 
tion. Since the thermally excited triplet states are separated from a non-magnetic ground 
state by an energy gap, the number of excitations and the magnetic susceptibility should 
decrease and tend to zero below the transition temperature. Both static susceptibility 
measurements and ESR-studies || |7|, § performed on crystals of pure CuGe0 3 showed 
that the susceptibility rapidly decreases below T$p =14.5 K. The transition temperature 
to the dimerised state was obtained from the initial decrease in susceptibility. 

CuGeOs is not a perfect one dimensional magnet, the intrachain exchange J c is larger 
but not much stronger than the interchain exchanges = 0.1 J c and J a = — 0.0 1J C 0. A 
next-nearest-neighbour exchange along the chains with a significant value of J' ~ 0.36J C 



also probably exists in this compound || [10], [U]]. If the crystal lattice of CuGe0 3 was 
harder, this crystal would be antiferromagnetic due to the intra- and interchain exchange 
interactions with a Neel temperature of about T N ~ (Jc-h) 1 ~ 10 K |T3|. But the spin- 



Peierls state wins the competition with long-range antiferromagnetic order. Doping by 
impurities makes long range antiferromagnetic order possible again. The presence of im- 
purities diminishes the temperature of the spin-Peierls transition and at lower temperatures 



a Neel ordering occurs |13|, [I4[ . The spin-Peierls dimerization and the antiferromagnetism 



were found to coexist in doped crystals. This phenomenon was explained by Fukuyama 



et al |L5[ and Khomskii et al |TTJ, who showed by different methods that the impurity or 
the defect in the dimerization should be surrounded by a region of antiferromagnetically 
correlated spins, forming magnetic clusters, or solitons. The tails of these solitons could 
overlap producing long range magnetic order. From such a point of view the reduction of 
the spin-Peierls temperature is proportional to the concentration of the dopant for small 
concentrations. 

The impurity or the point magnetic defect inserted into the nonmagnetic spin-Peierls 
matrix generates a "many spin" object consisting of the impurity itself and of several neigh- 
bouring Cu-ions spins. The magnetic object formed due to the presence of the impurity 
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is therefore a mesoscopic type object containing a number of spins which is intermediate 
between a microscopic and macroscopic systems, while the magnetisation remains micro- 
scopic. 

Single crystals of CuGe03 with the following substituting impurities were investigated: 
Si, Ti on Ge-sites 0, Zn, Mg 0, |T|], Ni and Co @ on the Cu-sites. 



The suppression of the spin-Peierls transition and occurrence of antiferromagnetic order 
at low temperatures are common features for different dopants. In the antiferromagnetic 
phase the magnetic moment per Cu-ion is strongly reduced. This reduction depends upon 
the type of dopant and their concentration. For the 3.2 % Ni-doped crystal the effective 
moment is /i e // = 0.16 ± 0.03/zb while it decreases to 0.06 ± 0.03/is for the 1.7% Ni-doped 
sample [pOj] . For the 3.2% Zn-doped crystal fj, e jf ~ 0.2fi B [p3| . The direction of the easy 
axis of the antiferromagnetically ordered state also depends on the dopant: the easy axis 
is directed along the c-axis for Zn 



Si [24], Co while for the Ni-doped crystals it is 



directed along the a-axis |2"T 



The aim of the present paper was the study of the ESR signals from the impurity seeded 
magnetic clusters and the search for the antiferromagnetic resonance (AFMR) in Ni-doped 
single crystals Cui-^Ni^GeOs. These were the same crystals which were investigated earlier 
by neutron scattering [20]. 



2 EXPERIMENTAL DETAILS AND SAMPLES 

In our experiments we used a set of ESR spectrometers with transmission type cavities 
in the frequency range 8-80 GHz. The microwave cavities were placed in a hermetically 
isolated volume immersed in a liquid helium bath and filled with a small amount of He-gas, 
enabling one to vary the temperature of the cavity containing the sample in the range 1.3 
- 20 K . The magnetic resonance absorption lines were recorded through the dependence 
of the transmitted microwave power on the applied magnetic field. The reduction of the 
transmitted power is proportional to the power absorbed by the sample when the absorbed 
power is low. In the paramagnetic state the intensity of the absorption integrated over 
magnetic field is proportional to the static susceptibility of the spin system. 

The magnetization curves were obtained using an Oxford Instruments vibrating sample 
magnetometer. 

Single crystals of Cui_ x Nij.Ge0 3 with x=0.017 and x=0.032 were produced by the 



crystal growth procedure described in |20| . These were the same samples used for neutron 
scattering experiments. These experiments showed the appearance of antiferromagnetic 
order at Tjv=2.3 K in the 1.7% Ni-doped crystal and at 4.2 K for x— 3.2% . The transition 
to the spin-Peierls state at Tsp = 11.5 K was observed for the 1.7% Ni-doped crystal, while 
for the sample containing 3.2%Ni the spin-Peierls transition was not clearly observed. The 
values of the exchange integrals Jf, = 0.7 ± 0.1 meV, J c = 1.8 ± 0.3 meV and of the spin 
wave energy gap A m 0.18 meV were obtained from the dispersion curves of the magnetic 
excitations in the 3.2% Ni-doped CuGe03 at T =1.5 K. To investigate the influence of 
the Ni-doping at low concentrations, when the dopant atoms do not interact, crystals of 
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Cui_ x Ni x Ge03 with x < 0.005 were grown by recrystallization of ceramic samples in air 
using a horizontal floating-zone method. The growth rate was 5-7 mm/h. The ceramic 
samples were prepared by annealing in air of nonstoichio metric mixture of CuO, GeC>2 and 
Ni 2 3 at 950° C for 24 hours. 



3 EXPERIMENTAL RESULTS 

3.1 The temperature evolution of the ESR line and the AFMR 
spectrum 

The evolution of the ESR line with temperature for the pure crystals of CuGeO"3 is described 
in || [7], The intensity of the resonance absorption diminishes rapidly below T$p. The 
nonmonotonic change of g-factors g a , gb and g c takes place at the temperature variation 
in the range between 14.5 and 4 K. This change occurs near the values which are close to 
g=2.1 and the magnitude of this change does not exceed 4 % for g a and g^ and 1 percent 
for g c . The additional line with g a = 1.82, g c =1.45, g^ =1.86 occurs below 4 K. The main 
line splits into three narrow lines at low temperatures. At a temperature of 4 K the ESR 
integral intensity of pure samples from our set of crystals is about 3xl0 -3 of the intensity 
of a paramagnet with one electron spin per Cu-ion. The residual ESR signal of the pure 
crystals is attributed to defects in the structure and to residual impurity content, as well 



as to the boundaries of the spin-Peierls domains p5 . 

The intensity of the ESR signals observed at low temperatures in the Ni-doped samples 
is larger than the intensity of the ESR in pure samples. The value of the intensity for the 
£=0.5% corresponds approximately to the concentration of the inserted impurities, while 
for x >1.7% the intensity is less than that of a paramagnet with the corresponding amount 
of S =1 spins. The evolution of the ESR line with temperature for samples with different 
concentrations of impurity is shown in Figures 1, 2, 3. 

In contrast to the pure material, the value of the g-factor strongly depends on the 
temperature below Tsp in Ni-doped crystals. The resonance field is shifted toward higher 
fields with decreasing temperature and g-factor diminishes by about 20%. For the smallest 
concentration, £=0.005, (Figure 1) an additional weak line with a g-factor of about 2 arises 
at low temperatures, indicating an ESR spectrum which is intermediate between the pure 
and doped samples. 

The sample doped with 1.7% Ni showed an analogous temperature shift of the ESR line 
(Figure 2). The linewidth has a maximum at a temperature of about 2.5 K (see Figure 4), 
which is close to the Neel temperature for this crystal. There is no observable shift of 
the line position between the Neel temperature 2.3 K and the lowest temperature of our 
experiments 1.3 K in the frequency range 18-60 GHz. 

The sample containing 3.2% Ni demonstrated a transformation from a single-line spec- 
trum into the spectrum typical of an orthorhombic antiferromagnet with several resonant 
lines. Three lines appear at this transformation when the magnetic field is parallel to the 
direction of spin ordering (Figure 3), and a gap opens in the absorption spectrum for other 
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directions of the magnetic field. This transformation occurs at 4 K corresponding to the 
Neel temperature for this sample. Above the Neel temperature, the resonance field slightly 
increases as the temperature decreases, but this shift is smaller than for more lightly doped 
samples. The temperature evolution of the resonance field is given in Figure 5. 

We performed measurements of the resonance field in the frequency range 18-75 GHz. 
The resonance frequency-field dependencies for 0.5% Ni and 1.7% Ni-doped samples were 
linear and gapless in the temperature interval between 1.3 and 20 K. A typical dependence 
of the resonance frequency / on the magnetic field H is given in Figure 6. The temperature 
dependencies of the g-factors taken at different frequencies are given in Figure 7. The 
sample doped with 3.2% Ni is presented here only at T > T/v, where the spectrum is 
gapless. The values of ^-factors are obtained by the relation: 

g = 2f/ 1 H with 7 = 28 GHz/T. 

Note that the deviation of the g-factors with temperature does not tend towards the 
value of the g-f actor for Ni 2+ ions which has the value about 2.3. At T > T$p there is 
also a remarkable deviation of g-factor for the doped crystals from the value of the pure 
material. This deviation increases with increasing concentration. 

The AFMR spectrum of the 3.2% Ni doped sample at T=1.8 K which has well pro- 
nounced gaps of 22 GHz and 33 GHz is presented in Figure 8. 



3.2 Spin-Peierls transition 

The drop of the integrated intensity of the ESR line marks the temperature of the spin- 
Peierls transition both for pure and doped crystals. The appropriate temperature depen- 
dence of the integrated intensity is shown in Figures 9-a,b. The transition temperatures 
obtained from the ESR data are given in the Table 1. 



Table 1 The temperatures of the spin-Peierls transition obtained from the ESR intensity vs 
temperature dependencies for the samples of Cu(i_ a; )Ni a ;Ge03 



X 





0.5% 


1.7% 


3.2% 


Tsp, K 


14.0±0.5 


12.5±0.5 


11.5±0.5 


8.0±1 



The values of T$p for doped crystals are smaller than those for pure CuGe03. The value 



of T S p for the 1.7% doped sample agrees well with the neutron scattering data pOR , while for 
the 3.2% doped sample the neutron scattering data cannot be used for the determination 
of Tsp because of the almost complete collapse of the spin-Peierls excitations. 



3.3 Magnetization curves 

The magnetization vs magnetic field was measured at different temperatures for the 1.7% 
Ni doped sample to investigate the magnetic properties of the sample in the temperature 
range of the strong ^-factor evolution. Figure 10 illustrates the M(H) dependencies at 
different temperatures for H \\ a. For the clear demonstration of the nonlinear contribution 
a fixed linear part is subtracted and the difference AM is plotted. Significant nonlinearity 
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arises at low (T<6 K) temperatures. The inset of this Figure demonstrates the spin-flop 
transition at H=0.3 T, H || a. For H || b and H \\ c the M(H) curves are similar except 
for the spin-flop steps, which are observable only for H || a. The nonlinear parts for the 
a and c-directions of the magnetic field are the same, for the 6-direction the nonlinear 
contribution is about 50% greater. 



4 DISCUSSION 

4.1 x-T phase diagram 

The values of the spin-Peierls transition temperature and the Neel points obtained from 
ESR intensity are in a good agreement with the x — T phase diagram for the Ni-doped 



samples presented in [26| with a linear decrease of T$p with x in the interval 1.7-2.9%. Our 
additional points at x =0.5% and 3.2% extend the experimental verification of the linear 
region. 



4.2 Antiferromagnetic resonance 

Doped crystals, which demonstrated the Neel order at low temperatures, show different 
types of magnetic resonance spectra. The sample doped with 1.7% Ni has the Neel point at 

2.3 K, but did not show a characteristic gap in the magnetic resonance spectrum (Figure 6). 
The sample doped with 3.2% Ni has a gap in the resonance spectrum (Figure 8). This 
spectrum may be described as the AFMR spectrum of an orthorhombic antiferromagnet, 
obtained by means of the Landau-Lifshitz equations of sublattice magnetizations motion 
within a molecular field approximation. Taking the axes a, c, b as easy, second easy and 
the hard axis correspondingly, we obtain the resonance frequencies v 12 as follows |27f| : 



H || a, H < Hgp \ 






K2/7) 2 = ^[(l + a 2 )# 2 + Ci + 


c 2 ± 




((l-a 2 ) 2 i/ 4 + 2(l + a) 2 (C 1 


+ C 2 )H 2 + (C 1 -C 2 ) 2 fi} 


(1) 


H\\ a, H > Hsf '■ 






K/7) 2 = H 2 — Ci 




(2) 


W7) 2 = C 2 -d 






H || c : 






(z/ 1 / 7 ) 2 = ^ 2 + C 1 




(3) 


(zy 2 / 7 ) 2 = C 2 






H || b : 






K/7) 2 = H 2 + C 2 




(4) 


(z/ 2 / 7 ) 2 = d 
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here C\ t 2 = 2H e H a i ja 2 (H e ,H a i jCl 2 are exchange field and fields of anisotropy respectively), 
Hsf = (ZHeHai/a) 1 / 2 is the spin-flop field and a = 1 — X\\/x±- The applied magnetic field 
H and the anisotropy fields are considered to be much smaller than the exchange field. 

By fitting our data according to equations (1-4) we obtained the following parameters 
of the AFMR spectrum: H SF =1.2±0.05 T, 7 = 24.6 GHz/T, 2H e H al = 0.85 ± 0.10 T 2 , 
2H e H a2 = 1.80 T 2 , a =0.75±0.2. The gap values for the second-easy and hard axis direc- 
tions of the magnetic field are 22 GHz and 33 GHz (0.11 meV and 0.17 meV) respectively, 
which are in agreement with the value of 0.18 meV obtained from neutron scattering |20] . 



The spin-flop magnetic field is characterized by the wide band of absorption instead 
of the single resonance frequency (see Figure 8). According to our data Hsf=^--2 T. This 
value is in good agreement with the value of the spin-flop field 1.1 T found in [21] from 
the magnetization curves for the Ni-doped crystal with x=0.033. 

From the exchange integral values Jb ~ 0.7 meV, J c m 1.8 meV described in Section 2 
and assuming S = 1/2 and H e = 4(J C + J b )S '/ '(s^b) we estimate the exchange field value 
H e 27 T. Further, from our data on the AFMR gap we deduce the value of the anisotropy 
fields: H al » 0.017 T and H a2 « 0.036 T. 

The single line ESR spectrum transforms into the AFMR spectrum with a gap for 
x =3.2% but for x =1.7% the transition to the Neel state is marked only by the maximum 
of the linewidth. If the conventional approach to the AFMR frequency derivation were 
valid for the sample containing 1.7% Ni, the observed value of the spin flop transition 
(Figure 10) would correspond to an AFMR gap of 7.5 GHz. For the frequency 18 GHz at 
H || a it should shift the resonance field by 0.05 T towards higher fields and for H \\ c to 
lower fields. This shift should be visible in our experiments by lowering temperature from 
T N =2.3 to 1.3 K. Nevertheless no shift exceeding 0.02 T was observed. This discrepancy 
between the static and dynamic properties of the impurity-induced antiferromagnetic state 
within the spin-Peierls matrix may be attributed to the low value of the order parameter 
for the Neel state. The sublattice magnetization equals only 0.06 of the nominal value. 
Usual procedure for deriving the AFMR frequency is therefore not valid because of the 
assumption that the spins are arranged in the form of hard sublattices, without a reservoir 
of disordered spins. In the case of the coexistence of the spin-Peierls and Neel ordered states 
the magnetically ordered part of the magnetization should interact with the disordered 
spin-Peierls background via the strong exchange interaction J c . It will probably cause an 
unusual type of magnetic resonance frequency or an overdamped mode. 

4.3 Spin clusters and magnetic resonance 

The value of g-factor in Ni-doped samples differs strongly from the corresponding values 
of pure crystals and of crystals doped with other impurities ^1, This difference takes 



place both above and below Tsp. Below the spin-Peierls transition there is also a strong 
anisotropy in the ^-factor. Because of the pronounced difference between the observed g- 
factor and that expected for individual Ni- and Cu-ions, we consider that clusters of several 
spins coupled by the exchange interaction are responsible for the discrepancy. The ^-factor 
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of the cluster of ions coupled by the symmetric Heisenberg exchange takes the averaged 



value between the ^-factors of isolated ions [28]. The observed g-factor value is well outside 



of this interval. The possible reason for this striking ^-factor deviation is the formation of 
clusters containing several spins coupled both by symmetric and antisymmetric exchange 



interactions [25]. 

The formation of clusters should occur in doped CuGeC>3-crystals around the dopant 
ions, as described in Sec. 1, because the defect is surrounded by several antiferro magnet- 
ically correlated spins. The characteristic length of the reduction of the correlated spin 



component on moving away from the impurity should be about 7 interionic distances |[T 



Clusters containing three S = 1/2 ions are known to display an effective ^-factor which 
is smaller than the ^-factor of an isolated ion. The anisotropy of the ^-factor for the cluster 
is larger than that of the isolated ions [BO, 31 , . This change in the ^-factor is described 
by taking into account the spin-orbital interaction combined with the Heisenberg exchange 
in the form of a Dzyaloshinski-Moriya antisymmetric exchange which is allowed when the 
symmetry of the pairs of the interacting ions is low enough. According to |29j the reduction 
of the effective ^-factor is of the order of D/5J. Here D is the antisymmetric exchange 
coefficient and S J is the difference between the exchange integrals within a triangular spin 
cluster. The presence of the Dzyaloshinski-Moriya term in pure CuGeOs was proposed 
[32| to explain the ESR linewidth at high temperatures. In addition, the Dzyaloshinski- 



m 



Moriya term describing the interaction of two neighbouring magnetic ions may arise below 
Tsp because of the lowering of the local symmetry resulting from the dimerization. 

Figure 11 shows an impurity atom embedded in the dimerised matrix. The dimerization 
is disturbed in the vicinity of the impurity. There are no symmetry centers for Ni-Cu pairs. 
The region around the impurity atom also does not have the symmetry center. Therefore 
Dzyaloshinski-Moriya interactions coefficients should have different nonzero values for each 
pair of ions in the vicinity of the impurity. 

We consider the models of the spin cluster with the final and relative small numbers 
of spins to describe qualitatively the ESR spectra and to evaluate the coefficients of the 
Dzyaloshinski-Moriya interaction. The larger the number of spins under consideration 
the more realistic is the model because the real cluster is formed on the basis of the 
impurity in the infinite chain. For the large number of spins in the model the correlated 
component of the spins lying far from the Ni-spin should be strongly reduced due to the 
dimerization. Therefore the magnetic properties of the model cluster should be independent 
of the number of spins. The possible values of the total spin of the cluster arising due to 
the substitution of one Cu-ion per Ni-ion are S=l and 5=1/2. In order to obtain S =1 
the spin of the Ni-ion must be uncompensated and therefore the Cu-ions are divided into 
dimers in another way when compared to the undisturbed chain. The correlation of the 
dimerization in the neighbouring chains will be violated in this case. For the case of the 
total spin S =1/2 the spin of the Cu-Ni pair is uncompensated. This total spin value does 
not cause the rearrangement of the dimers lattice except for the single pair Cu-Cu replaced 
by Cu-Ni pair. Therefore the total spin S =1/2 of the cluster corresponds to the lower 
energy of the perturbation of the dimerised lattice and the 5* =1/2 models with the finite 
number of spins are probably more realistic then the models with S = 1. 
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We have analyzed the five- and six-spin- models for the cluster. The five-spin- model has 
S =1 in the ground state and the six spin model has S =1/2. The three-spin cluster would 
be the simplest case, but it has S = in the ground state and is therefore nonmagnetic. 
The four spin model could not include dimers from both sides of the Ni-ion. 

The five- and six-spin-clusters are shown schematically in the Figure 11 and are de- 
scribed by the Hamiltonian 

K K i=K 

^ 9iaf^BH a Si a (5) 
i,j i,j i=l,a=a,b,c 

Here gi a are g-tensor components of Cu and Ni-ions, fig is the Bohr magneton. are 
the nearest-neighbour and the next-nearest-neighbour exchange integrals and Dy are the 
vectors of the antisymmetric exchange. The different pairs of ions are taken only once in 
the sums ■ K is the number of spins equal to 5 or 6. 

In the following analysis we shall assume the Djj vectors to be parallel to each other and 
perpendicular to c axes of the crystal. For the five spin cluster the ions are labeled CU5- 
CU4-M1-CU2-CU3 and the nonzero Jjy are taken as follows. J54 = J23=10.6 meV are the 
Cu-Cu exchange integrals; J 12 = Ju =5 meV are the exchange integrals between Ni and Cu 
ions, they are estimated as a half of the Cu-Cu exchange integral [j33fl ; J42=3.6 meV is the 
next-nearest-neighbour exchange in the pure material ||); J i3 and J 15 , the next-nearest- 
neighbour exchange of the Ni-ion are taken to be 2 meV. The vectors are supposed to 
be parallel to each other and perpendicular to the chains. 

The diagonalization of the Hamiltonian (|5]) gives the energy levels for S=l states of the 
five-spin-cluster depending on the Dij values and magnetic field. The levels characterized 
with 5*2 = 0, ±1 are spilt by the terms of (|5]). The transitions with the momentum 
change of ±h are therefore separated by a gap of 1.0 meV in zero field. The formulae 
determining the gap and the values of g-factor are given in the Appendix. This gap is the 
main feature differentiating the S=l and 5=1/2 cases of five- and six-spin models. 

The lowest S =1/2 states of the of the six-ion-clusters are also separated by a gap 
which depends on the exchange integrals and Dy values. A magnetic field splits these 
levels, this splitting depends on magnetic field and D^. However, the transitions between 
the sublevels with the ±h momentum change remain gapless. The modification of the ESR 
spectrum at low frequencies is restricted here to the renormalization of the effective g- 
factors. The energy levels are obtained by the diagonalization of the energy matrix 10x10 
for the S =1/2 states of the six-spin-cluster. For the cluster with the ions labelled as 
Cu 5 -Cu 6 -Ni 1 -Cu2-Cu 3 -Cu4 with the following values of the exchange integrals Ju=5 meV, 
Ji 6 =5 meV, J 2 3=9.8 meV, J 34 = J 56 =10.6 meV the terms dominating in g-factors at 
magnetic field perpendicular and parallel to Dzyaloshinski-Moriya vectors are: 



4gm - 9cu 4g Ni + hg Cu (4Z?i 6 + 3D, 



9± * 3^ 9^ 54EI (6) 



-21 



9\\ » 5 (7) 
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where 



E2I ~ J56 ~ 77^61 + 77^62 — 77^15 (8) 
3 6 9 

is the energy interval between the ground state and the lower excited state of the cluster. 

To obtain g±=1.6 as observed in our experiments we have to assume Diq m ^3meV. 
The exact diagonalization of the matrix gives the close values of \D{j\ = 3.4meV. 

The values of the averaged spin projection at the site for six-spin-cluster obtained by 
the described procedure are given in the Table 2 and are shown schematically in Figure 11. 
The following numerical values are used at this procedure in addition to the exchange 
integrals given above: D 12 = D 16 =D23=D 3 4=D 65 =3.4meV. 

Table 2. Average values of the spin projections of the ions constructing a model six-spin 
cluster 





s 5 


s 6 


Si 


s 2 


s 3 


s. 




see text 


0.16 
0.117 


-0.118 
-0.125 


0.614 
0.591 


-0.145 
-0.152 


0.007 
0.013 


-0.017 
-0.02 



Thus the results obtained on the basis of the six-spin-model correspond to the observed 
gapless ESR spectrum with the anisotropic g-factor deviating from the free spin g-factor. 
The five-spin-model does not correspond to the observed ESR signals due to the absence 
of a gap in the observed spectrum. Therefore the spectra obtained confirm the proposed 
S =1/2 structure of the cluster. 

Note that in our experiments we observed the reduction of g-factor with temperature 
for all principal orientations of the magnetic field, while the model with collinear vectors 
predicts the deviation for only one principal direction, the ^-factors for the other two 
principal directions should be close to 2. Nevertheless these simple models demonstrate 
the possible mechanism for the reduction of the value of effective ^-factor and of its strong 
anisotropy. The deviation from the free spin ^-factor value of the other components of the 
^-factors tensor may be probably provided by any noncollinearity of the vectors Djj. 

The nonlinearity of magnetization curves at low temperatures which is shown in Fig- 
ure 10 also confirms the formation of spin clusters with their intrinsic degrees of freedom, 
because the energy levels of cluster are separated by the gaps depending in a nonlinear way 
on magnetic field. The existence of a nonlinear susceptibility was reported as evidence for 
Dzyaloshinski-Moriya interactions in three- ion clusters in |Q. The growth with temper- 
ature of the linear part of the susceptibility, visible in Figure 10, is obviously due to the 
temperature dependence of the concentration of the triplet excitations of the spin-Peierls 
state. The nonlinear contribution to the magnetization curves diminishes and vanishes at 
the temperature rise from 1.6 to 6 K. At the temperature about 9 K the magnetization 
curve becomes again nonlinear but with the convexity directed down. The differential sus- 
ceptibility grows with magnetic field. This behaviour may be explained by the destruction 
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of the spin-Peierls state by applied magnetic field p5[. This destructive magnetic field 



could be of the moderate value when the temperature is close to Tsp- 

The nonlinear magnetization curve at 2 K might be described as a sum of the linear 
M(H) curve and of the magnetization of a paramagnet with the concentration of 1CT 3 
S —1/2 ions per one Cu-ion, the nonlinearity being ascribed to the paramagnetic satura- 
tion. However the analogous curve for 6 K could not be described in this way because the 
characteristic field of the saturation is not shifted to higher fields. 

The change of the g-factor of Ni-doped samples with temperature may be ascribed to 
the freezing out of the reservoir of the triplet excitations of the spin-Peierls state below 
T,sp- It results in the switching-off of the exchange narrowing and in the breaking of the 
collective character of the precession mode of the impurities and triplet excitations of the 
spin-Peierls state. This process is similar to the temperature evolution of the resonance 



spectrum of the triplet excitations in ion radical salts with a singlet ground state [36 



At low temperatures only the impurity mode (spin-cluster mode) survives, which leads to 
the unusual value of the g-factor. At intermediate temperatures the intermediate ^-factor 
is observable, due to the exchange mixing of the spins states. 



5 CONCLUSION 

The data obtained correspond to the x-T phase diagram show the linear dependence of 
Tsp and T^r with doping concentration. From the AFMR spectra we found that the im- 
purity induced antiferromagnetism in Cuo.96sNio.o32Ge03 can be described in terms of a 
molecular field theory as a conventional orthorhombic antiferromagnet with the easy axis 
directed along a, hard axis along b and second-easy axis along c. From AFMR spectra we 
obtained gaps of 22 and 33 GHz and anisotropy fields (H a i ~ 0.017 T, H a2 ~0.036 T).The 
antiferromagnetic ordering of the 1.7% Ni-doped sample could not be considered within the 
molecular field approximation because of the small value of the sublattice magnetization 
accompanied by the magnetically disordered spin-Peierls background. This fact results in 
the absence of a conventional AFMR spectrum. The anomalous value and temperature 
dependence of g-factor in the spin-Peierls state of Ni-doped crystals indicates the formation 
of spin clusters around the doping magnetic ions. The presence of magnetic clusters with 
their own internal degrees of freedoms is confirmed by observation of nonlinear magneti- 
zation curves in the magnetically disordered state. Further theoretical and experimental 
investigations of the cluster structure are necessary for a more detailed interpretation of 
ESR spectra and magnetization curves. 
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7 Appendix 

The calculations of the energy levels of five ion linear cluster Cu 5 -Cu4-Ni 1 -Cu 2 -Cu3 are 
based on the Hamiltonian || The energies of three low laying triplets with S=l are given 
approximately by 



where 



= ^ + ^ T y (It - I 2 a - ^) 2 + 2(J U + / 13 ) 2 - (9) 
E 3 = ^(h + I 2 )-I 23 . (10) 



h = -\{J& + Jab)- (11) 

h = -~(Jl4 + Jis)- ^23 + ^45- (12) 
/ 3 = ~^(Jl2 + ^13) - ^45 + ^23 (13) 

I12 = ~^2^ U ~ ^ 15 )' = "^Tf^ 12 ~~ ^ 13 ^' ^ 23 = 4^ 24 ' 

The Dzyaloshinski-Moriya interaction splits the ground state triplet into the levels with 
the energies E s , E a ^ 



E s = Ei, E a fi = Ei — C 2 - — ^ 14 ^ ■ (15) 



E 3 — Ei 

where 



G12 = ^(D 12 + D 23 -D 13 ), G14 = ^(D 14 + D 45 -D 15 ). (16) 

r= 0.5(J 2 + / 3 ) + / 2 3-^i 

1 {[0.5(/ 2 + J 3 ) + / 2 3- J Ei] 2 + 0.5(/ 12 + / 1 3) 2 } 1 /2- ^ > 
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The energies in the magnetic field are determined by the secular equation 



E s - e g x eff (3H x 
g x eff /3H x E a -e 
gl ff (3H z -iglffpH, 








(18) 



Here g x ff , g y eff ,and g z ef} are given by 




where 



-(/i2 + /i 3 )/v / 2 



(21) 



{[0.5(/ 2 + J 3 ) + J 23 - £i] 2 + 0.5(/i 2 + /i 3 ) 2 } 1/2 ' 



The Dzyaloshinsky-Moriya vector G12 — G14 is supposed to be parallel to the y-axis. 
When the external magnetic field is perpendicular to the vector G12 — G14 and E s — E a » 
gf3H the frequency-field dependence is quadratic. 

References 

[1] Hase M, Terasaki I, Uchinokura K 1993 Phys. Rev. Lett. 70 3651 

[2] Pytte E 1974 Phys. Rev. B 10 4637 

[3] Nishi M, Fujita O, Akimitsu J 1994 Phys. Rev. Lett. 50 6508 

[4] Regnault L P, Ain M, Hennion B, Dhalenne G, Revcolevschi A 1996 Phys. Rev. B 53 



[5] Hirota K, Cox D E, Lorenzo J E, Shirane G, Tranquada J M, Hase M, Uchinokura K, 
Kojima H, Shibuya Y, Tanaka I, 1994 Phys. Rev. Lett. 73 736 

[6] Oscroff S, Cheong S-W, Fondado A, Aktas B, Fisk Z 1994 J. Appl. Phys. 75 6819 

[7] Honda M, Shibata T, Kindo K, Sygai S, Takeuchi T, Hori H 1996 J. Phys. Soc. Jpn 
65 691 

[8] Smirnov A I, Glazkov V N, Vasil'ev A N, Leonyuk L I, Coad S M, Paul D McK, 
Dhalenne G and Revcolevschi A 1996 JETP Lett. 64 305 



5579 



13 



Riera J, Dobry A 1995 Phys. Rev. B 51 16098 

Braden M, Wilkendorf G, Lorenzana J, Ain M, Mclntyre G J, Behruzi M, Heger G, 
Dhalenne G, Revcolevschi A 1996 Phys. Rev. B 54 1105 

Khomskii D, Geertsma W, Mostovoy M 1996 Proceedings of the 21-st International 
Conference on Low Temperature Physics, Prague 1996 Chech. J. Phys. 46 Suppl. S6 
3239 

Hennessy M J, McElwee C D, Richards P M 1973 Phys. Rev. B 7 930 

Regnault L P, Renard J P, Dhalenne G, Revcolevschi A 1995 Europhys. Lett. 32 579 

Lussier J-G, Coad S M, McMorrow D F Paul D Mck 1995 J. Phys.: Cond. Matter 7 
L325 

Fukuyama H, Tanimoto T, Saito M 1996 J. Phys. Soc. Jpn. 65 1182 

Weiden M, Richter W et al 1996 Physica B 225 177 

Hase M, Uchinokura K et al 1996 J. Phys. Soc. Jpn. 65 1392 

Hase M, Sasago Y, Uchinocura K, Kido G, Hamamoto T 1995 JMMM 140-144 1691 

Hase M, Hagiwara M, Katsumata K 1996 Phys. Rev. B 54 R3722 

Coad S, Petrenko O, Paul D McK, Lussier J-G, McMorrow D F 1997 Physica B 239 
350 

Koide N, Sasago Y, Masuda T, Uchinokura K 1996 Proceedings of the 21-st Inter- 
national Conference on Low Temperature Physics, Prague 1996. Czech J. Phys. 46 
Suppl.S4 1981 

Anderson P E, Liu J Z, Shelton R N 1997 Phys. Rev. B 56 11014 

Sasago Y, Koide N, Uchinokura K, Martin M C, Hase M, Hirota K, Shirane G 1996 
Phys. Rev. B 54 R6835 

Nojiri H, Hamamoto T, Wang Z J, Mitsudo S, Motokawa M, Kimura S, Ohta H, 
Ogiwara A, Fujita O, Akimitsu J 1997 J. Phys. Condens. Matter 9 1331 

Smirnov A I, Glazkov V N, Leonyuk L I, Vetkin A G, Eremina R M 1998 submitted 
to JETP 

Coad S, Lussier J-G, McMorrow D F, Paul D Mck 1996 J. Phys. Condens. Matter 8 
6251 

Nagamiya T, Yosida K, Kubo R 1955 Advances in Physics 4 Nu.13 



14 



Abragam A, Bleaney B. Electron Paramagnetic Resonance of Transition Ions, Claren- 
don Press, Oxford 1970, Chapter 9. 

Belinskii M I, Tsukerblat B S, Ablov A V 1974 Fix. Tverd. Tela (Sov.Phys. Solid State) 
16 989 

Yablokov Yu V 1978 J. Mol. Struct. 46 285 

Yablokov Yu V, Voronkova V K, Mosina L V Paramagnetic resonance of exchange 
magnetic clusters (in russian) Moscow, Nauka, 1988 

Yamada I, Nishi M, Akimitsu J 1996 J. Phys. Cond. Matter 8 2625 

Eremin M V, Rakitin Yu V 1977 Physica status solidi 80 579 

Tsukerblat B S, Novotvortsev V M, Kuyavskaya B Ya, Belinskii M I, Ablov A V, 
Bazhan A N, Kalinnikov V T 1974 Pisma Zh. Exp.Teor.Fiz (Sov.Phys. JETP Lett ). 
19 525 

Hase M, Terasaki I, Uchinokura K, Tokunaga M, Miura N, Obara H 1993 Phys. Rev. 
B 48 9616 

Chesnut D B, Phillips W D 1961 J. Chem. Phys. 35 1002 

Conncll H M Mc, Griffith H O, Pooley D 1962 J. Chem. Phys. 36 2518 



15 



Figure captions 

Figure 1. Temperature evolution of the ESR line for the 0.5% Ni-doped sample. H || c, 
/=36.7 GHz. The arrow marks the resonance field of the free electron spin (g=2). 

Figure 2. Temperature evolution of the ESR line for the 1.7% Ni-doped sample. H || b, 
/=36.0 GHz. The arrow marks the resonance field for g=2. Narrow line is DPPH-mark. 

Figure 3. Temperature evolution of the ESR line for the 3.2% Ni-doped sample. H || a 
(easy axis direction). The arrow marks the resonance field for g=2. 

Figure 4. Temperature dependence of the ESR linewidth for the 1.7% Ni-doped sample, 
/=36.4 GHz. 

Figure 5. Temperature dependencies of resonance fields for the 3.2% Ni-doped sample. 
/ ~36 GHz. The arrows mark the resonance fields for g=2. 

Figure 6. The ESR spectrum of the crystal doped with 1.7%-Ni H || c, T=1.8 K 
Figure 7. Temperature dependencies of ^-factors. 

Figure 8. Spectrum of the antiferromagnetic resonance of Cu .968Nio.o32Ge0 3 at T—1.8 K. 

Figure 9-a,b. Temperature dependencies of the integrated ESR intensity for pure and 
doped crystals. Dotted curve is the the susceptibility of a paramagnet containing 0.5% of 
spins S=l, g =2.3. 

Figure 10. Magnetization curves for the 1.7% Ni-doped sample at H || a. 

Figure 11. The structure of the dimerised lattice around the impurity atom. The five- 
and six-spin-clusters. The arrows indicate the direction and the value of the average spin 
projections within the six-spin-cluster with Dy =0. 
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